January 14, 2013 0:59 WSPC/INSTRUCTION FILE 

cleymans' worku ' mpla ' v2 



The Hagedorn Temperature Revisited. 



J. Cleymans 

UCT-CERN Research Centre and Department of Physics, University of Cape Town, 
Rondebosch 7701, South Africa 

D. Worku 

UCT-CERN Research Centre and Department of Physics, University of Cape Town, 
Rondebosch 7701, South Africa 

The Hagedorn temperature, Th is determined from the number of hadronic resonances 
including all mesons and baryons. This leads to a stable result Tg = 174 MeV consistent 
with the critical and the chemical freeze-out temperatures at zero chemical potential. We 
use this result to calculate the speed of sound and other thermodynamic quantities in the 
resonance hadron gas model for a wide range of baryon chemical potentials following the 
chemical freeze-out curve. We compare some of our results to those obtained previously 
in other papers. 
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1. Introduction 

In 1965 Hagedorn ^ proposed that the number of hadronic resonances increases 
exponentially with the mass m of the resonances. The idea, which was debated 
strongly when first proposed , has since been widely accepted and discussed by many 
authors I 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |1T] xhe concept was based on the assumption that the 
observed increase in the number of hadronic resonances would con tinu e towards 
higher and higher masses as more experimental data became available ^^21, xhe scale 
of the exponen tial increase deter mines the value of the Hagedorn temperature, Th- 



Group to revisit the original analysis of Hagedorn to update the value of Th- 
This resulted in a surprising wide spread of possible values, with large variations 
as to whether one considers mesons or baryons with values ranging from Th = 141 
MeV to Th — 340 MeV depending on the parametrization used and on the set of 
hadrons (mesons or baryons). There thus exists uncertainty as to the value of the 
Hagedorn temperature. These have two origins: 

• sparse information about hadronic resonances certainly above 3 GeV, 

• the analytical form of the Hagedorn spectrum, especially the factor multi- 
plying the exponential. 



Recent 
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The first item will probably never be resolved satisfactorily due to the width of 
resonances and also due to their large number making it difficult to identify them. 
Splitting the spectrum into baryons and mesons further decreases the quality of the 
fits to the mass dependence of the mass spectrum. We therefore propose to stick to 
the original analysis of Hagedorn and consider a sum over all resonances, baryons, 
mesons, strange, non-strange, charm, bottom etc.. This is the state that is produced 
at the Large Hadron Collider (LHC), namely, a hadronic ensemble containing all 
possible resonances. The result is shown in Fig. [1] and leads to a good determina- 
tion of Tu because the range in m is reasonably large extending up to 3 GeV before 
reaching a plateau presumably due to the parsity of hadronic resonances above this 
value. Details about the parametrization used will be presented below. 
The Hagedorn temperature naturally leads to the notion that hadronic matter can- 
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Fig. 1. Cumulative number of hadronic resonances as a function of m. Again the hadronic data 
are made up of all resonances, including baryons, mesons and also heavy resonances made up of 
charm and bottom quarks. 

not exceed a limiting temperature and increasing the beam energy in proton— proton 
and proton— antiproton collisions results in more and more hadronic resonances 
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being produced without a corresponding increase in the temperature of the final 
(freeze-out) state. This is the situation observed at the highest energies at the 
LHC. It was suggested long ago that the Hagedorn temperature is connected to 
the existence of a different phase in which quarks and gluons are no longer confined. 
At present, the Hagedorn temperature is often understood as the temperature of 
the phase transiti on from hadrons to a quark-gluon plasma. 

A recent analysis l ^^ l ^^ l investigated the hadron resonance gas model to 
show the stability of various thermodynamic quantities in heavy-ion collisions when 
the Hagedorn spectrum is used literally, i.e. without a cut-off on the number of 
resonances beyond a certain mass. In particular, the authors explored the addi- 
tion of the hadron resonance gas model including an exponentially large number 
of undiscovered resonances which are naturally included in the Hagedorn model. 
Their results showed that the hadron resonance gas model gave different results for 
thermodynamic quantities but the overall chemical analysis was reasonably stable. 

Furthermore, the use of so-called Hagedorn States (HS), based on the expo- 
nentially increasing spectrum, close to t he critical temperature can explain fast 
chemical equilibration by HS regeneration 1221 and provide a unique method to com- 
pare lattice res ults for Tc using thermal fits HS provide a lower than thermal 
fits without HS These authors estimated effects by extending the hadron mass 
spectrum beyond 3 GeV for Th = 200 MeV and assume that high mass excited K* 
states p rod uce one kaon while producing multiple pions thus further reducing the 

In this paper, we extend previous work^^SESEZI j^y extending explicit expressions 
of relevant thermodynamic quantities for non-zero chemical potentials. In particular 
the speed of sound, , can be considered as a sensitive indicator of the critical 
behavior in strongly interacting matter. The results show a sharp dip of in 
the critical region, which is an indication that thermodynamics in the vicinity of 
confinement is indeed driven by the higher excited hadronic states. 
The outline of this paper is as follows. In section 2, we present the influence of the 
Hagedorn spectrum on the hadron yields to find the thermodynamic parameters and 
explain the basic concepts used in this paper. In section 3, we derive the number, 
energy and entropy densities and also the speed of sound. In section 4, we show 
the results using the hadron and its extension resonance gas model for particular 
thermodynamic quantities and discuss the relationship of with temperature and 
chemical potentials and compare thermodynamic quantities, like energy and entropy 
density. The last section covers the conclusion. 



2. Motivation 

The particle data table contains hundreds of hadronic resonances including the 
well known stable hadrons like nucleons, pions, hyperons (w, S, S), kaons etc.. In 
our calculations we used a list of hadronic resonances including in total 250 baryons 
and antibaryons, 348 mesons and antimesons (counted without considering their 
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isospin and spin degeneracies). Most of the hadronic resonances decay quickly via 
strong interactions before reaching the detector, hence they are usuaUy identified 
via their decay products. The mass of a decaying particle is equal to the total energy 
of the products measured in its rest frame. 

The basic idea of this paper is to add resonances using the Hagedorn model for the 
spectrum. Using the hadronic data we show the relationship between the number 
of hadronic resonances and the mass in Fig. [T] where we took hadrons with masses 
up to 11.019 GeV. The density of states obtained this way has been fitted using the 
following equation ' ■ ; 

Ph{m) = ^ ^ ^ . ./4 exp 1;^ ), (1) 
where c and ttiq are constant parameters given in the table below 



Parameter 




c (GeV^)3/2 
mo (GeV) 
Th {GeV) 


0.16 ± 0.02 
0.5 

0.174 ± 0.011 



3. The Hadron Resonance Gas Model and its Extension 

The thermodynamic properties of the Hadron Resonance Gas Model (HRGM) can 
be determined from the partition function 

In Z{V,T,fi) — J dm [p M (m) In Zb{m,V,T, fi) 

+ pB{m)\nZf{m,V,T,fj.)], (2) 

where the gas is contained in a volume V, has a temperature T and chemical 
potential /x, Zb and Zf are the partition functions for an ideal gas of bosons and 
fermions respectively with mass m, pMim) and ps(m) are the spectral density of 
mesons and baryons. By using Eq.[21 one can compute the number denisty n, energy 
density e, entropy density s, pressure P, speed of sound and specific heat Cy. 
hadron properties enter these models through pm,b- The HRGM model takes the 
observed spectrum of hadrons up to some cutoff of mass A, defined by 

< A 

PALsim) = ^ gi5{m-m,), (3) 

i 

where mi are the masses of the known hadrons and gi the degeneracy factor. In order 
to explore the stability of results obtained using the HRGM, variant of these models 
is often used in which one takes the observed spectrum of states up to a certain cutoff 
of mass A and above this one includes an exponentially rising cumulative density 
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of hadron states, which is defined in Eq. [TJ This defines the Extended Hadron 
Resonance Gas Model (EHRGM) . The density of states becomes 

Ph{m) = V gib{m - rrii) + — — exp ( — ) e{m - A), (4) 

where the model parameter c and are fitted to data on the cumulative distribu- 
tion of the sets of hadrons, h. Typically this model uses mo = 0.5 GcV. Moreover, 
the parameters are determined from the hadronic spectrum with masses up to 3 
GeV as shown in Fig. [TJ The results for c and Th are given in Table 1. 



4. Derivation of thermodynamic quantities in EHRGM 

We are considering here the Boltzmann distribution for simplicity in order to present 
some of the results in a compact form. The partition function for a single particle 
is given by 

where K2 is modified Bessel function. The meson mass distribution is taken to be 
given by 



mi<A 



1=1 



m 

(m2 + ^2)5/4 



4.1. Derivation for the particle densities and pressure 

The particle density can be written as the sum of the two terms 

n ^ HM +nB, (7) 

where hm — ^t^Zm/V and ns = InZs/V are number densities of mesons and 
baryons respectively. They are defined by 

„ mi<\ 
1 _ ( P 

nM,B 



^ ( Pi\ \ M,B 2 / "^A I 

1=1 



(8) 



where Hi = Sifis + BifiB + QiPQ, for our case we consider an isospin symmetric 
system, where /xq — 0.0 GeV. The pressure is given by 

P-T^+T^=P,,+Ps, (9) 

where Pm — Tum and Pb = Tub are the pressure of mesons and baryons respec- 
tively. 



I 



I 
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4.2. Derivation for energy density 

The energy density is given by 

T''d\nZM,B 

£m^B - -y of ' y^^) 

where Ki is the modified Bessel function, em and £b are energy densities of mesons 
and baryons respectively. 



where 



Ah^c 



m \ m 
2T^e^P ^ jK2[7f )dm. 



A 



H 



(11) 



From the above expressions, we can obtain the entropy density 

e + P - ns^is - nsiiB 
f ' 

where ns and ng are the net number densities for strange and baryonic particles 
respectively 

_ T d\nZ 

ns{B) — Tj-^ • 

Vdfis(B) 



(12) 



(13) 



4.3. Derivation of the speed of sound 

In hydrodynamic models the speed of sound plays an important role in the evo- 
lution of a syst em and is a n ingredient in the understanding of the effects of a 
phase transition | 25|26|27 |28j jg well-known that the speed of sound has to be 
calculated at constant entropy per particle (s/n). This makes the calculation more 
complicated than the one at zero chemical potential where it is s ufficient to keep 
the temperature fixed. In our extension we take the condition into account for 
non-zero baryon and strangeness chemical potentials, imposing overall strangeness 
zero. The squared speed is thus calculated starting from 

The complete expression for the speed of sound can be rewritten as 



(14) 



Cf(T,/.) 



'dP} 
KdT) 



) + (C) (^) + (c) (^) 



(15) 



where the derivatives and can be evaluated from two conditions. The first 



condition comes from the requirement that the ratio (s/n) must be kept fixed. 



- = o, 

n 



(16) 



I 



I 



January 14, 2013 0:59 WSPC/INSTRUCTION FILE 

cleymans' worku ' mpla ' v2 



The Hagedorn Temperature Revisited. 7 

while the second condition fohows from the conservation of strangeness 

ns = rig. (17) 
The resulting expressions for and are presented in detail in the appendix. 

5. Results using HRGM and EHRGM 

There is not much difference between the HRGM and EHRGM models at low tem- 
perature, T <C Th since the heavy resonances do not play an important role there. 
Alternatively, at high temp eratures, T ^ Th, it should agree with the results of the 
lattice simulations of QCD l ^^ l ^^ l ^^ l. In our case, we have determined the temper- 
ature and chemical potential dependencies of the pressur e aiid the number, energy 
and entropy densities using the chemical freeze-out curve ^l ^^l^^ i Using these re- 
sults, we have then calculated the speed of sound. The thermodynamic variables 
obtained using HRGM, we plot the energy and entropy densities scaled by the ap- 
propriate powers of T, which is shown in Fig. [2^, [3^ and H^, at the Th = 0.174 
GeV didn't observe any sudden change of the thermodynamic variables e/T^ and 
s /T^ close to the critical temperature, it showed smooth shape as the temperature 
goes beyong Th. 

Moreover, when we use the EHRGM, we see a different behaviour at the critical 
temperature which is shown in Fig. ^jp, [Sb and IHd. These show that if the system 
undergoes a first-order phase transition, both temperature and pressure remain con- 
stant as hadronic matter is converted from hadron gas to quark-gluon plasma; the 
energy and entropy density, however, change discontinuously. This leads to show us 
a new state of matter. 

The value of the speed of sound remains well below the ideal-gas limit for mass- 
less particles = 1/3, even at very high temperatures, energy densities and pres- 
sure. In Fig. 1313 represents, when it crosses the transition point in the case of full 
QCD, we expect as before that should vanish beside that is inversely pro- 
portional to the specific heat Cy, and this can lead to a divergence at the critical 
temperature Tc- Of course, due to finite volume effects, the velocity of sound will 
most likely not completely vanishes at T = Tc- The velocity of sound versus the 
temperature as shown in Fig. 0)3 is to evaluate dP/de, we have followed!^ and first 
expressed P and e in physical units, using EHRGM in Eq. 01 The temperature and 
chemical potential dependence of all the relevant thermodynamic quantities shows 
unique behaviour at the critical point Tc, specially for the speed of sound, there is 
a pronounced dip as evidence for the phase transition in the system. Based on our 
calculations, the velocity of sound at T = Tc is different from zero; it is expected to 
become zero but due to finite volume effects. 

6. Conclusions 

In this paper we have made a new analysis of the numbe r of hadronic resonances 
using the latest information from the Particle Data Group This leads to a tem- 
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Fig. 2. (a) This figure show that the energy density e in units of T* calculated on the HRGM as 
a function of the temperature T. (b) This figure show that the energy density e in units of T** for 
EHRGM as a function of the T. The full— lines are the result of the EHRGM that accounts for all 
mesonic and baryonic resonances. 



perature which is consistent with the most r ecent r esults based on lattice QCD 
estimates of the phase transition temp erature 1^ 01^ and also the chemical freeze- 
out temperature at zero baryon density I 19 | 20 | 21 | have extended calculations of 
the speed of sound to non-zero baryon and strangeness chemical potentials keeping 
s/n fixed I^Sl. This is done for both the hadronic resonance gas model (HRGM) and 
the extended hadronic resonance gas model (EHRGM) which includes the exponen- 
tially increasing spectrum of hadrons following the Hagedorn parametrization. The 
EHRGM shows that the speed of sound goes to zero at the phase transition point 
while the HRGM shows a smooth dip followed by an increase. 
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Fig. 3. (a)This figure show that the entropy s in units of T'^ calculated in the HRGM as a function 
of the T. (b) Entropy s in units of T^ using EHRGM as a function of the T. The full— lines are 
the result of the EHRGM that accounts for all mesonic and baryonic resonances. 



Appendix A. Appendix: Speed of sound at non-zero chemical 
potentials. 



The speed of sound is given by ^ 



dP 



(A.l) 



where s/n is the entropy per particle which is kept fixed, using the variables T,jiB 
and this can be rewritten as 

liP.\ -L ( dP_\ ( di£s A 1 ( J)P_\ ( d4iB_\ 

CUT,,)^ ^'"^' y'''l)''^( (A.2) 



de \ I dp.B 
dfiB J \ dT 



where the derivatives and can be evaluated using two conditions. The first 
condition comes from keeping the ratio {s/n) constant. From the derivative one 
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Fig. 4. (a) Showing the speed of sound (T, jj,) calculated on the HRGM as a function of the T 
and /i. (b) The speed of sound C|(T, n) in EHRGM with resonance mass truncation, m < 2 GeV 



obtains 
which implies 

nds = sdn. 

In terms of T, /is and fis this equation can be written as 
n (^) dT + n (^) d„s + n (|^) d^. 



(A.3) 



(A.4) 



dn 
df 



dn 



= s I t::^ ] dT + s [ ) dfiB + s - — d^is- 



dm 



dn 



\dlJLS 



(A.5) 
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divide the above expression by dT on both sides, so that it becomes 
ds\ f ds \ f d^iB 

.9tJ+"'— , ' — 

'diiB\ , ^( dn\ f dpLs \ 



dfiB 
dn 
df 



V dT 
f dn 



( ( ^ 



dT 



+ s 



dfis J \dT J 



(A.6) 



Rearranging the above expression in order to write ^^j^ in terms of one obtains 



ds 



df) 



dn 
d^iB 

df) 



/ dnB \ 
[dT ) 
ds 
dfJ.s 



f dn 



( 

\dT ) 



Defining 



A^n\ — \-s\ — 



B = 



ds_ 

df 
ds 
d^B 
ds 
dus 

The final expression for condition one is 

d^B 1 



dn 
df 
dn 



C 



dfiB 
dn 



dT 



B 



A + C 



dT 



The second condition comes from overall strangeness neutrality, which is 

ns = ns 



(A.7) 

(A.8) 
(A.9) 
(A.IO) 

(A.ll) 
(A.12) 



where ns and ng are the strange and antistrange particle densities. Similarly the 
derivative of equation Eq. IA.12I should thus satisfy 



d{ns) = d{ns) 

this implies equation Eq. IA.13I can be expressed as 

' dns \ / d^iB \ f dns 
dtiB)\dT) \d^ls 



(A.13) 



/ dns^ 
V dT 
dr 



dT 



dni 
dfiE 



dfiB 

~fr 



( dr 



dT 
dT 



(A.14) 



We can apply the same method as above to write 
relation 

'dng 



in terms of ^ for the above 



/ dns 
\dfj.B 
dng 
dT 



dfiB 
'dT 

fdns\ 
Ws) 



dfis 



dT 



(A.15) 
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we define L — ns, hence it represents that the number of strangeness density for 
baryons and mesons 

L = n^ + n¥. 



and R — rig, the number of antistrangeness density for baryons and mesons 

i? = n| + nf 



Define now 



dT dT 

Hence, the final expression for condition two becomes 



^ dL OR 

E^^ 7^, (A.16) 

(A.17) 

D^'^~'^. (A.18) 



dT 



(A.19) 

(A.20) 
(A.21) 



dT E 

Finally, by equating equation lA.llI and IA.19I wc find 

d^is _AE-BD 
It ~ BF -CE' 

and 

diiB CD- AF 
~dT ~ BE -CE' 
Which is the relation used in the text. 
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